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ABSTRACT 

The stability of compressible two- and three-dimensional boundary layers is 
reviewed. The stability of 2D compressible flows differs from that of incompressible 
flows in two important features: There is more than one mode of instability 

contributing to the growth of disturbances in supersonic laminar boundary layers and 
the most unstable first-mode wave is three-dimensional. Whereas viscosity has a 
destabilizing effect on incompressible flows, it is stabilizing for high supersonic 
Mach numbers. Whereas cooling stabilizes first-mode waves, it destabilizes second- 
nic^e waves. However, second-order waves can be stabilized by suction and favorable 
pressure gradients. The influence of the nonparallelism on the spatial growth rate 
of disturbances is evaluated. The growth rate depends on the flow variable as well 
as the distance from the body. Floquet theory is used to investigate the subharmonic 
secondary instability. 

1. INTRODUCTION 

The aim of this paper is to review the state of the art of the stability of 
compressible boundary layers. The study discusses the influence of Mach number, 

Reynolds number, cooling, suction, pressure gradients, wave angle, and 
nonparallelism. Subharmonic secondary instability is also discussed. 

The earliest attempt at formulating a compressible stability theory was made by 
Kuchemann (ref. 1) who neglected viscosity, the mean temperature gradient, and the 
curvature of the mean velocity profile. Lees and Lin (ref. 2) and Lees (ref. 3) were 
the first to derive the basic equations for the linear parallel stability analysis of 
compressible boundary layers. This theory was extended by Dunn and Lin (ref. 4), 

Reshotko (ref. 5), and Lees and Reshotko (ref. 6). These early theories were 
asymptotic or approximate in nature and proved to be valid only up to low supersonic 
Mach numbers. The use of direct computer solutions to exploit the full compressible 
stability equations was initiated by Brown (ref. 7) and Mack (ref. 8). An extensive 

treatment of the parallel stability theory for compressible flows is given by Mack 

(refs. 9-15). As the Mach number increases, the dissipation terms become important 
and a three-dimensional disturbance cannot be treated by an equivalent two- 
dimensional method as is usually done for the incompressible case. Mack (refs. 
10,15) found that neglecting the dissipation terms can lead to a 10% error in the 
disturbance amplification rate. 

It is an interesting facet of compressible two-dimensional boundary layers that 
the most unstable first-mode wave need not be parallel to the freestream as the Mach 
number approaches one. At supersonic speeds the most unstable first-mode wave is 
oblique or three-dimensional. 

The most important feature of the stability of supersonic laminar boundary 
layers is that there can be more than one mode of instability contributing to the 
growth of the disturbance. The first mode is similar to the Tol lmien-Schl ichting 
instability mode of incompressible flows, while the second and higher unstable modes 
are unique to compressible flows. Mack (ref. 10) found that there are multiple 
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values of wave numbers for a single disturbance phase velocity whenever there is a 
region of supersonic mean flow relative to the disturbance phase velocity. For 
incompressible flows, higher modes are associated with higher wave numbers at 
different phase speeds. In contrast with the first mode, the most unstable second 
mode is two-dimensional. As the Mach number increases to the hypersonic regime, the 
second mode displays growth rates that are higher than those of the three-dimensional 
first mode. However, the maximum growth rate is less than that of the first mode at 
zero Mach number. 

These stability theories treat the mean flows as quasiparallel flows. Some 
incomplete attempts to account for the nonparallel flow effects by including either 
the normal velocity or some of the streamwise derivatives of the mean flow were given 
by Brown (ref. 16), Gunness (ref. 17), and Boehman (ref. 18). Complete nonparallel 
theories for two-dimensional flows were developed by El-Hady and Nayfeh (ref. 19) and 
Gaponov (20) and for three-dimensional flows by Nayfeh (21) and El-Hady (22). The 
growth rate in a parallel flow is independent of the flow variable and the distance 
from the wall, whereas the growth rate in a nonparallel flow (growing boundary layer) 
depends on the flow variable and the distance from the wall. This complicates the 
interpretation of experimental data for comparison with the results of stability 
theory. 

In contrast with the case of incompressible flows, rigorous stability 
experiments are very difficult at supersonic speeds because (a) the spatial and 
temporal resolution of instruments at supersonic speeds is less than those at low 
speeds, (b) with the exception of Kendall, experimenters have less control and 
knowledge of the disturbances at supersonic speeds, and (c) the interference of 
traversing probes at high speeds is due to the high aerodynamic loads which 
necessitate strong and bulky walls. Therefore, most of the information on the 
stability of high-speed flows is macroscopic rather than microscopic. The term 
macroscopic refers to measurements of the onset of turbulence and the extent of the 
transition region, whereas the term microscopic refers to measurements of the 
evolution in space and time of the fluctuations present in the flow that are 
sufficient for the identification of the instabilities that lead to transition and 
the validation of the proposed theoretical models. It should be noted that 
macroscopic experiments are difficult to relate directly to stability theory, whereas 
microscopic experiments, which provide information about the unstable disturbances 
and their growth, can be better related to stability theory, which studies the 
development of individual components of the disturbances corresponding to a certain 
frequency or a wavepacket. Thus, microscopic experiments that use controlled 
disturbances are more desirable for corroboration with theory than experiments that 
study natural disturbances arising from one source or another in the boundary 
layer. The natural disturbances represent a set of space and time components, 
whereas controlled experiments can provide disturbances with a given frequency and a 
given spanwise wavenumber. 

Whereas experimenters developed various credible techniques to introduce 
controlled artificial disturbances in incompressible boundary layers, the technique 
of Kendall seems to be the only credible technique at supersonic speeds. Laufer and 
Vrebalovich (ref. 23) used a high-speed valve and Demetriades used a siren mechanism 
attached to a flat plate to introduce their artificial disturbances. The valve 
opened and closed a narrow slit in the surface of the plate to allow periodic air 
pulses of certain frequency to disturb the boundary layer. Kendall introduced small 
artificial disturbances by a glow discharge between two electrodes embedded in the 
surface of the flat plate skewed at a specified angle to the spanwise direction. 
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thereby introducing disturbances with a specified wave angle. The rest of the 
available experiments were performed using natural disturbances. 

Almost all measurements reported on the experimental stability of boundary 
layers were made by means of either hot-wire anemometers or hot films following 
disturbances in these boundary layers. The hot-wire response is a combination of 
velocity, density, and temperature fluctuations. The hot-wire response when operated 
at high constant overheat is proportional to the mean square of the mass-flow 
fluctuations. To describe the disturbance fully, one needs measurements of all 
fluctuation characteristics such as the rms amplitude, spectra, and propagation speed 
as functions of both y and x. Almost all reported experiments for compressible flows 
measured |pu| at various x-stations by placing hot wires or hot films at a transverse 
location where the mean-flow conditions are the same. Laufer and Vrebalovich (ref. 
23) reported measurements at different constant y/L positions, while Kendall (refs. 
24,25), Demetriades (refs. 26-28), Lebiga et al (ref. 29), and Stetson, Thompson, 
Donaldson and Siler (refs. 30-33) reported measurements at one constant y/L located 
in the wideband energy peak. 

Laufer and Vrebalovich (ref. 23) carried out measurements of the neutral 
stability curves, amplification rates, wavelengths, and amplitude distributions at 
the Mach numbers 1.6 and 2.2. They performed their measurements in the JPL 20" 
supersonic wind tunnel where the freestream turbulence level was reduced to about 1% 
by means of damping screens. Laufer and Vrebalovich performed their measurements for 
natural as well as artificial disturbances. The stability characteristics of natural 
disturbances in supersonic flows at Mach numbers between 1.6 and 8.5 were examined by 
Kendall (refs. 24,25) in the JPL 20" supersonic wind tunnel. In some of these 
experiments, the side walls of the tunnel were turbulent and hence radiated sound. 
Mack (ref. 12) tried to compare the free oscillations in the parallel stability 
theory with Kendall's data. The comparison was satisfactory only for the case M = 
4.5. In an attempt to account for the response of the boundary layer to the incoming 
sound waves, Mack (ref. 12) included a forcing term at the first neutral stability 
point and found a better agreement with Kendall's data. The same characteristic 
features of the boundary- layer response to the incoming sound waves were observed by 
Lebiga et al (ref. 29) in their experiments at M^ = 2. Oemetriades (refs. 26,27) 
presented experimental results for hypersonic bounci ary- layer flows. He studied the 
streamwise amplitude variation of both natural disturbances and disturbances 
artificially excited with a siren mechanism attached to a flat plate. 

2. PROBLEM FORMULATION 

In this paper, we consider the linear quasi-parallel as well as the nonparallel 
stability of two- and three-dimensional compressible boundary layers. Moreover, we 
consider the linear secondary instability of two-dimensional primary waves in a two- 
dimensional compressible boundary layer. The basic equations for the linear 
stability analysis of parallel-flow compressible boundary layers were first derived 
by Lees and Lin (ref. 2), Lees (ref. 3), and Dunn and Lin (ref. 4), using the small 
disturbance theory. For excellent references on the compressible parallel stability 
theory, we refer the reader to the papers of Mack. For the nonparallel theory of 
compressible boundary layers, we refer the reader to the papers by El-Hady and Nayfeh 
(ref. 19) and Nayfeh (ref. 21). 

Lengths, velocities, and time are made dimensionless using a suitable reference 
length L*, the freestream velocity U*, and L*/U*, respectively. The pressure is made 
dimensionless using p*U* 2 , where p*”is the freestream density. The temperature, 
density, specific heats,” viscosity^ and thermal conductivity are made dimensionless 
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using their corresponding freestream values. The gas is assumed to be perfect. 
Since the pressure is constant across the boundary layer 


p 


T 

m m 


= 1 and p„ 
m 



( 2 . 1 ) 


where the subscript m refers to mean flow quantities, y is the ratio of the specific 
heats of the gas, and M^ is the freestream Mach number. 

To formulate the problem so that the disturbance equations can be specialized to 
all these cases, we assume the basic flow to be a time-dependent three-dimensional 
flow and superpose on it small disturbances to obtain total flow quantities of the 
form 


q(x,y,z,t) = q b (x,y,z,t) + q(x,y,z,t) 


( 2 . 2 ) 


Here q b (x,y,z,t) stands for a basic state quantity and q(x,y t z,t) stands for a small 
unsteady disturbance. Substituting the total flow quantities u, v, w, p, p, y, and t 
into the Navier-Stokes equations, subtracting the basic state, and linearizing, we 
find to first order that the disturbance equations are given by 
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The constants r and m are given by 


r = ^ (e + 2) and m = -| (e - 1) 


( 2 . 10 ) 


where e = 0 corresponds to the Stokes hypothesis. The Reynolds number R and Prandtl 
number Pr are given by 
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(2.11a) 


U*L* 


R = 




Pr = 


P*C 


where v*, u*, and <* are the freestream kinematic viscosity, dynamic viscosity, and 
thermal” conductivity, respectively, and ci is the dimensional specific heat at 
constant pressure. We assume that the dynamic viscosity is a function of the 
temperature only so that 



The boundary conditions at the wall are 


(2.11b) 


u = v = w = T = 0 at y = 0 (2.12) 

The boundary conditions on the velocity fluctuations u and w represent the no-ship 
conditions and the boundary condition on the velocity fluctuation v represents the 
no-penetration condition. For a gas flowing over a solid wall, the temperature 
remains at its mean value unless the frequency is small (i.e., stationary or near 
stationary disturbances). The boundary conditions as y - ® are 


u(y), v(y), w(y) , p(y), and T(y) are bounded as y - » (2.13) 

As will be described later, neutral subsonic disturbances decay to zero 
as y - ®, whereas neutral supersonic disturbances do not vanish as y - ®. 

3. QUASIPARALLEL PRIMARY INSTABILITY 

In this section, we consider the three-dimensional stability of a steady three- 
dimensional boundary layer. In general, the mean flow in a boundary layer varies 
with the streamwise coordinate x and spanwise coordinate z. However, at high 
Reynolds numbers, this variation is small over distances the order of the wavelengths 
of the disturbances. Consequently, most stability analyses neglect the streamwise 
and spanwise variations of the mean flow, the so-called parallel-flow assumption. 
Thus, the basic flow is approximated by 


“b - V y) ’ v b " °> w b * T b ' V y >- p b ‘ P m (x) f 3 ' 1 * 

Using the parallel-flow assumption reduces equations (2. 3)-(2.9) to a system of 
linear partial differential equations whose coefficients vary only with y. 
Consequently, the variables t, x, and z can be separated using the so-called normal 
mode assumption 


[u,v,p,T,w] = [^(y), c 3 (y), c 4 (y), c 5 (y), c,(y)| exp[i(ax + bz - wt)] (3.2) 

where a and 6 are the wave numbers in the streamwise and spanwise directions, 
respectively, and u> is the frequency. For a temporal stability, a and 6 are real 
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but ui is complex. For spatial stability, o> is real but a and g are complex, 
general case, a, 8, and w are complex. 

Substituting equations (3.1) and (3.2) into equations (2.3)-(2.9), (2 
(2.13) yields the eigenvalue problem 
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Introducing the Stuart transformation (ref. 34) 


kn j = aC; + 8 c 7 

(3.10) 

kn 7 = ac 7 - 

(3.11) 


where k = /a 2 + e 2 does not reduce the eighth-order system (3.3)-(3.9) into a sixth- 
order system due to the coupling in the energy equation due to a single dissipation 
terms, which couples the energy equation to the other equations for 3D waves in 2D or 
3D boundary layers. Neglecting this term reduces the eigenvalue problem from an 
eight-order complex problem to a sixth-order complex problem, thereby resulting in a 
large saving in computations. The error introduced to this simplification is a 
function of the Mach number, Reynolds number, and wave angle. For the case of 
insulated flat plates. Mack found that the maximum error is less than 5 X for all Mach 
numbers at R = 1500. However, most stability calculations being formed are based on 
the eighth-order rather than the sixth-order system. 

Two methods of solution have been employed (refs. 11,35,36). The first reduces 
the system of equations (3.3)-(3.7) into an eighth-order system of ordinary 
differential equations, determines an exact solution outside the boundary layer that 
satisfies the boundary conditions (3.9), uses the resulting linearly independent 
solutions as initial conditions to integrate the first-order equations to the wall, 
employs a Gram-Schmidt orthonormalization scheme to keep the solutions linearly 
independent, and uses an iterative scheme such as Newton-Raphson procedure to 
determine the desired complex eigenvalue. The second method uses a finite-difference 
scheme to reduce the governing equations into a system of algebraic equations that is 
solved using standard techniques to determine the desired eigenvalue. 

According to the first approach, the eigenvalue problem (3.3)-(3.9) is converted 
into a system of first-order equations by letting 


5 2 = Du = Dcj, c 6 = DT = Dc 5 , and c 8 = Dw = Ds 7 (3.12) 
Then, the eigenvalue problem becomes 

Dc = F(yK (3.13a) 
?! = z 3 = = c 7 = 0 at y = 0 (3.13b) 
5 are bounded as y - ® (3.13c) 


where ^ = {c { ,c 2 , . . . 8 } and the elements a^ .• of the matrix F are given in Appendix 
A. For a given R and a mean flow U m , W m , and T m , the eigenvalue problem (3.13) 
provides a dispersion relation of the form 


to = to(a,6) 


(3.14) 
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In general, w, a, and 6 are complex. Thus, equation (3.14) provides two relations 
among the six real parameters cy, o^, e r , 6-j, u> r , and When four of these 
parameters are specified, equation (3.14) can be solved for the remaining two 
parameters. For the case of temporal stability, a and 8 are assumed to be real and 
fixed, and equation (3.14) provides w = w + iw.. Thus, it follows from equation 
(3.2) that the temporal growth rate is w.. r For ihe case of spatial stability, u> is 
assumed to be real and fixed, and equation (3.14) provides two relations among the 
remaining four parameters a r , c^, e r , and 8-j. The parameters a and e define a real 
wavenumber vector £ whose magnitude k is given by r 


2 2 U 

k = (cy + B r ) k 


and whose direction (wave angle) <i* is given by 


(3.15a) 


= arctan (e r /a r ) (3.15b) 

Moreover, the parameters a. and 6- define a growth vector a whose magnitude o is 
given by 1 1 

o = (a 2 . + e 2 .) H (3.15c) 

and whose direction (growth direction) ii> is given by 


= arctan (e^/a..) (3. 15d) 

In general, the growth and wave directions need not coincide and two more relations 
need to be specified to complete the formulation of the spatial stability problem. 
Using the method of multiple scales (refs. 37,38), Nayfeh (refs. 21,39) found that 
the ratio of the group velocity components must be real, thereby providing a third 
relation. Maximizing the total growth rate can provide the fourth needed relation. 
These points are discussed further in Section 4. Nayfeh and Padhye (ref. 40) used 
the complex group velocity vector to relate the problems of temporal and spatial 
stabilities. 

In the case of 2D spatial stability, the two additional relations are given by 
6=6+ ie. = 0 and Eq. (3.14) can be used to determine the complex parameter a for 
any given w. For the case of 3D spatial stability in boundary layers that depend on 
x and y only (such as 2D mean flows and flows over an infinite span swept wing), the 
analysis in Section 4 shows that 6 = constant and hence Eq. (3.14) can be used to 
determine the complex parameter a for any given w and e. In this case, the growth 
direction is often taken in the x-direction (i.e., 8j is assumed to be zero). 

To account for the slow growth of the boundary layer, one can improve the 
parallel-flow assumption by calculating the local stability of the mean-flow 
profiles, thereby obtaining values for a and 8 that vary slowly with x and z. Thus, 
equation (3.2) is replaced with 


u,v,p,T,wj = |c,(y), c 3 (y), c 4 (y), c 5 (y), c 7 (y) ]exp(ie) 


(3.16) 
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where 


(3.17) 


30 / x 30 / N 30 
3X = a ( X ’ Z )’ = 6 ( x ’ z )» it = ■ “ 

and the z, are governed by the eigenvalue problem consisting of equations (3.3)- 
(3.9). Tnis is the so-called quasi-parallel assumption. As discussed earlier, the 
growth rate vector is given by 


a = a.i + 6 . j 

It is independent of the transverse direction and the flow variable whose growth is 
investigated; neither of these statements is valid for the case of growing boundary 
layers as discussed in Section 4. 

3.1 Inviscid Case 

In contrast with the case of incompressible boundary layers, compressible 
boundary layers, even on flat plates, have inviscid instability, which increases with 
increasing Mach number. The inviscid instability is governed by equations (3.3)- 
(3.9) with R being set equal to infinity; that is, it is governed by 
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Equations 

(3.18)- 
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can be 

combined into 

a second-order 

equation 


governing z , 3 or c 4 . To accomplish this, we use the transformation (3.10) and 
(3.11). Thus, we add a times equations (3.19) to 8 times equation (3.21) and obtain 

ik(aU + bW - u>)n, + z, ,D(oll + eWj + ik'T r = 0 (3.25) 

' m m ' 1 3 v m nr m * 

Using equation (3.22) to eliminate s 5 from equation (3.18) and using equation (3.10), 
we obtain 
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ikrij + Oq 
Eliminating r\ l 
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M 


(3.28) 


where 


aU +SW -to 

m m 


(3.29) 


and 


M = 
r 


(aU+eW -0))M 
m m ° 

k/T 


(3.30) 


can be interpreted as the local Mach number of the mean flow in the direction of the 
wavenumber vector 1< = at + 8J relative to the phase velocity w/k. In general M r is 
complex and it is only real for neutral disturbances. In terms of x and M r , equation 
(3.20) can be rewritten as 


2 2 

ik M 

0^4 = “ 2 X 

M 


(3.31) 


Eliminating from equations (3.28) and (3.31) yields the following second-order 
equation governing x : 


M 


D x + D[ln — ^y]D x - k (1 - M ) x = 0 
1-M r 

Eliminating x from equations (3.28) and (3.31) yields 


(3.32) 


D 2 c 4 - D(lnM 2 )D? 4 - k 2 (l - m'h 4 = 0 


(3.33) 


In the freestream, T_, = 1 , U_ = 1 and W_ = 0 are constants and hence M r = M^ is 
a constant and equations (3.32) ana (3.33) reduce to 


D 2 x - k 2 (l - M 2 ) x = 0 (3.34) 

and 

D 2 £ 4 - k 2 (l - M 2 )s 4 = 0 (3.35) 

It is clear from equations (3 ? 34) and (3.35) that neutral disturbances decay in the 
freestream if and only if Mf < 1; these disturbances are termed subsonic waves. 
When Mf > 1, neutral disturbances do not vanish in the freestream and they represent 
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sound waves or Mach waves of the relative flow^ they are termed supersonic waves and 
may be outgoing or incoming waves. When Mf = 1, disturbances are termed sonic 
waves. The above classification is due to Lees and Lin (ref. 2). 

For the case of two-dimensional inviscid waves in a two-dimensional boundary 
layer, Lees and Lin (ref. 2) established a number of conclusions for the case of 
temporal waves. Their main conclusions are 

(i) The necessary and sufficient condition for the existence of a neutral subsonic 
wave is the presence of a generalized inflection point y > y in the boundary layer 
at which 


D(p DU ) = 0 (3.36) 

' m nr 

where 

W - 1 - r < 3 - 37 > 

CD 

The phase velocity c $ = ui/a s of this neutral wave is U m (y s ), the mean velocity at the 
generalized inflection point y $ , which is larger than 1 - 

oo 

(ii) A sufficient condition for the existence of an unstable subsonic wave is the 
presence of a generalized inflection point at y s > y 0 ; its phase velocity 

c > 1 - 1/M^. Compressible boundary- layer flows over insulated flat plates always 
have inflection points and hence they are unstable to inviscid disturbances. 

(iii) There is a neutral subsonic wave having the wavenumber a = 0 and the phase 
velocity c = c Q . 

(iv) When M* < 1 everywhere in the boundary layer, there is a unique 
wavenumber ct s corresponding to the phase velocity c $ for the neutral subsonic wave. 

Using extensive numerical calculations. Mack (refs. 11,15) established the 
existence of an infinite sequence of discrete wave numbers a 5p , corresponding to an 
infinite sequence of discrete modes when M r > 1 somewhere in tne boundary layer. He 
referred to the modes that are additional to the mode found by Lees and Lin as higher 
modes. In contrast with the first mode whose existence depend^ on the presence of a 
generalized inflection point, the higher modes exist whenever M r > 1, irrespective of 
the presence or absence of a generalized inflection point. The lowest Mach number at 
which the higher modes exist in the boundary layer on an insulated flat plate is 
2.2. It turns out that this is also the lowest Mach number at which subsonic higher- 
mode disturbances exist. The lowest of the subsonic modes is called the second mode 
and it is the most amplified of the higher modes. 

Later, Mack developed a simple theory that provides an approximation to the 
infinite sequence of wave numbers a . He neglected the second term in (3.32) and 
expressed the solution of the resulting equation that satisfies the boundary 
condition v(0) as 


x ~ - sin [ k sn f /M £ - 1 <Jy], y < y a 

0 

and the solution that decays as y - °° as 


(3.38) 
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(3.39) 


x * - 1 exp[- k $n J /I - M2 dy], y > y a 

Ya 2 

where y a is the value of y where M =1. He chose the constant in (3.39) to be -i so 
that the pressure is real and positive for y > y . He argued that D x must go to zero 
as y + y a because the pressure is continuous and finite at y = y a , thereby obtaining 
the eigenvalues 

k sn ^ / fV r ~ I dy = (n - A)*, n = 1, 2, 3, ... (3.40) 

0 

We note that the second term in equation (3.32) is singular at y a and hence the 
expansions (3.38) and (3.39) are invalid at y = y a . To connect these expansions, one 
needs to consider the neighborhood of y = y a , determine an expansion there, and match 
it to both (3.38) and (3.39). Alternatively, the Langer transformation can be used 
to determine a single uniformly valid expansion. Even such a consistent expansion is 
valid only when M r does not vanish anywhere, (i.e., there is no critical layer) and 
hence there is no generalized inflection point; otherwise, one cannot neglect the 
second terms in (3.32) and (3.33). An expansion valid when M- vanishes somewhere is 
determined below using the Olver transformation. Mack (refs. 11,15) called the waves 
when = 0 inflectional neutral waves and those which occur when M * 0 non- 
inflectional neutral waves. His numerical results for the case of non-inf lectional 
neutral waves (Figs. 1 and 2) are qualitatively in agreement with his simple theory, 
whereas his calculations for the case of inflectional waves (Figs. 3 and 4) differ 
from those obtained from his simple theory. According to his theory, s 4 (s) is 
positive for all modes, there are no zeros in the interval y > y a , and the number of 
zeros in y < y a increases by one for each successive mode. 

To determine a single uniformly valid expansion, we eliminate the second term in 
equation (3.33) using the transformation 


= M* 


(3.41) 


and obtain 


2 K' 2 K 2 

■ 1 ' r . r 


- (k (1 - M r ) - jjp— + ]♦ = 0 


(3.42) 


For the case of non-inf lectional waves (i.e., M * 0 everywhere), we use the Langer 
transformation (refs. 37,38) and obtain the approximate solution 


A 


1/4 


4/1— M2 

r 


c x A i ( c ) + c 2 Bi(0 


(3.43) 


and hence 


M r $ 


r 1 / 4 


V1-M2 

r 


CjAiU) + c 2 Bi(0 


(3.44) 


where Ai and Bi are the Airy functions of the first and second kind, respectively, 
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(3.45) 


| 5 , 3/2 » k / /I - N* dy 

y a 

and y = y, is the turning point at which M (y ) = 1. Since 

ex P [-f e /z \ 

Ai(t) ' Z — TTt as c ♦ ® (3.46) 

2/tt £ 


Bi (0 


exp(| 

A 


the vanishing of the disturbance as y - °° demands that c 2 
condition v( 0 ) = 0 or c 3 ( 0 ) = 0 demands that 


(3.47) 


Then, the boundary 


Ai' = 0 at y = 0 

(3.48) 

But 

Al (0 ‘ _ i /4 sink (-0 / + 4 it) as j- - - - 

/it (-0 ' 

(3.49) 

Hence 


y i 

cos[ k J VM2 - ldy+^7i] = 0 
0 

(3.50) 

Consequently, 


k sn J 3 /M2 - 1 dy = (n - fb, n = 1, 2, 3, ... 
0 

(3.51) 

whose righthand side differs from that of Mack by n. Moreover, the eigenfunction 
(3.44) with c ? = 0 differs from that of Mack in the presence of the factor 
M r /Vl - M 2 and the phase ^ u in the interval y < y a . 

For the case of inflectional waves (i.e., M r = 
expansion (3.43) is not valid as y -* y c and we need an 
neighborhood of y = y c . Using the Olver transformation 

0 at y = y c , y c > yj, the 
expansion that is valicfin the 
(ref. 37), we find that 

U y /l-M? dy)"* y 

♦ - C {b 1 I 3/2 [k J /I - M 2 dy] 

V1-M2 ' y c 


+ b 2 K 3/3 [k / /I - Hf dy)) 

X- 

(3.52) 
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which is valid in the interval y > y a . Since 

I 3/ 2 (0 ^ -= and K 3/2 ft) - /Jr as 5 - » (3.53) 

1 /2n ' 

the boundary condition that c 4 tends to zero as y ® demands that b x = 0. This 
expansion needs to be matched with the expansion (3.43) that is valid at the turning 
point y a . Expanding (3.43) for large k and for y > y a , expanding (3.52) for large k 
and for y < y c , means that bj = 0, we obtain 

_ k/y C /FMTdy c -kf* /ftM^dy c -kj* /T^dy 

/Jr b 2 e = ^_ e a + — e a (3.54) 

c 2/u ✓w 


Hence, c 2 = 0 and 

c/k y 

b 2 - — exp(- J c A - dy] (3.55) 

n/2 y a 

Imposing the boundary condition at the wall that v(0) = 0 or c 3 (0) = 0 yields the 
eigenvalues as in (3.51). Thus, the difference between inflectional and non- 
inflectional waves is in the shape of the eigenfunction. 

Using extensive numerical computations, Mack established a number of conclusions 
regarding supersonic stability: 

(i) In contrast with incompressible stability theory, there is more than one mode of 
instability and it is one of the additional modes that is the most unstable. Figure 
5 shows the variation of the maximum temporal amplification rate of 2D waves with 
Mach number. It is clear that below M = 2.2, the boundary layer on an insulated 
flat plate is virtually stable to invTscid 2D waves and that above M^ = 2.2 the 
second mode is the most unstable mode. Moreover, the maximum amplification rate 
increases sharply as M^ increases beyond 2.2. Furthermore, above M^ = 5, the first 
mode is not even the second most unstable mode. 

(ii) In contrast with incompressible stability theory, 3D first modes are more 
unstable than their corresponding 2D waves. However, 3D second modes are more stable 
than their corresponding 2D waves. Figure 6 shows the variation of the temporal 
amplification rate of the first and second modes with frequency for M^ = 4.5 and 
several wave angles. It clearly shows that the most unstable first-mode “wave has a 
wave angle that is approximately 60° and an amplification rate that is approximately 
twice the maximum amplification rate of its corresponding 2D wave. 

(iii) Whereas cooling can stabilize first-mode waves in accordance with the pre- 
diction of Lees (ref. 3), cooling destablizes second-mode waves. Figure 7 shows the 
variation of the maximum temporal growth rate of the first mode at M = 3.0, 4.5, 

and 5.8 and the second mode at M = 5.8 with the ratio of the wall temperature T w to 
the recovery temperature T f . ft is clear that 2D and 3D first-mode waves can be 
completely stabilized by cooling, whereas second-mode waves are destablished by 
cooling. 
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3.2 Viscous Stability for Flat Plates 


Again, using extensive numerical calculations. Mack (refs. 10,11) investigated 
the influence of Mach number on the viscous instability of supersonic flows past flat 
plates. He found that viscosity is stabilizing for both 20 and 3D first-mode waves 
when M > 3.0 and for second-mode waves at all Mach numbers; that is, the maximum 
ampl ifTcation rate (over all frequencies, and wave angles in case of 3D waves, at 
constant Reynolds number) decreases with decreasing Reynolds number. This result was 
disputed by Wazzan, Taghavi, and Keltner (ref. 41) who did not find a transition from 
viscous to inviscid instability with increasing Mach number but found that viscous 
instability persists to M = 6.0. Mack (ref. 42) reconfirmed his calculations for 

the case of temporal stability and obtained spatial stability results that agree with 
his earlier conclusions on the influence of viscosity on compressible stability. 
Moreover, the spatial stability calculations of El-Hady and Nayfeh (ref. 19) and 

Nayfeh and Harper (ref. 43) agree with those of Mack for at least three significant 
figures for all the calculations they performed. Moreover, the calculations of Malik 
qualitatively agree with those of Mack. 

Figure 8 shows the neutral stability curves of 2D wav^s calculated by Mack at 
M^ = 1.6, 2.2, 2.6, 3.0, and 3.8. He plotted a vs. R” to emphasize the higher 

Rlynolds number region. Comparing the neutral curve for M^ = 1.6 with that of the 

Blasius flow (i.e., = 0) shows that, although they have” the same general type, 

compressiblity drastically reduced the viscous instability, resulting in much lower 
neutral wave numbers. As the Mach number increases beyond 1.6, the viscous 

instability continues to weaken although the effect of the increase in inviscid 

instability continues to extend to lower and lower Reynolds numbers. When M^ reaches 
3.8, the viscous instability disappears and viscosity acts only to damp” out the 
inviscid instability. 

As in the case of inviscid instability, the most unstable first-mode waves are 
3D waves. Figure 9 shows the variation of the maximum temporal amplification rate of 
2D and 3D waves with Reynolds number for M = 1.3, 1.6, 2.2, and 3.0 and Figure 10 
shows the variation of the maximum spatiaf amplification rate of 2D and 30 waves 
at ip = 50° with R for M^ = 1.6 (ref. 19). The most unstable wave angles of 3D waves 
are shown in Figure 9. At M = 1.3, the instability of 2D and 3D waves is due to 
viscosity. As M increases, Ihe viscous instability decreases for both 2D and 3D 
waves. Whereas increasing the Mach number results in a drastic decrease in the 2D 
maximum growth rates, it produces only a slight change in those of the 30 waves. 

As in the inviscid case, the numerical results of Mack suggest that the 2D 
second- and higher-mode waves are more unstable than their corresponding 3D waves. 
Moreover, the maximum growth rate of second-mode waves drops sharply as the wave 
angle increases from zero. 

The lowest Mach number at which Mack was able to calculate 2D second-mode waves 
is M = 3.0 at which the minimum critical Reynolds number is 13,900. As pointed out 
in tie preceding section, the inviscid instability increases rapidly with increasing 
Mach number and hence one would expect the minimum critical Reynolds number to 
decrease rapidly to lower values as the Mach number decreases. In fact, Mack found 
that the minimum critical Reynolds numbers drops to 235 as the Mach number increases 
to 4.5. Moreover, at high Mach numbers second-mode waves have much higher growth 
rates than oblique first-mode waves. 

Whereas cooling can stabilize 2D and 3D first-mode waves, second-mode waves 
cannot be stabilized by cooling. On the contrary second-mode waves are destablized 
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by cooling. Malik (ref. 43) studied the influence of cooling on oblique first-mode 
waves at = 60° for M ro = 2 and 4.5 and second-mode waves for M ot = 4.5; in these 
calculations R = 1500 and the stagnation temperature is 560°R. Figure 11 shows that 
the oblique first-mode wave at M = 2 is completely stabilized when T w /T . = 0.7, 
whereas that at M^ = 4.5 is stabfized only when T /T, d = 0.48. On the other hand. 
Figs. 11 and 12 show that the second-mode wave is aestablished by cooling. In fact, 
the maximum growth rate increases rapidly with cooling. Malik found that the 
frequency of most amplified first-mode wave decreases with cooling whereas that of 
the most amplified second-mode wave increases with cooling. 

Malik (ref. 43) also investigated the influence of favorable pressure gradients 
(Fig. 13) and suction (Fig. 14) on the stabilization of second-mode waves. He found 
that each of them shifts the band of unstable frequencies to higher values and 
reduces the peak amplification. Consequently, it appears that, whereas cooling 
cannot be used to stabilize second-mode waves, they can be stablized using either 
suction or wall shaping to produce a favorable pressure gradient. 

Preliminary results of Malik (ref. 43) indicate that real gas effects tend to 
destabilize hypersonic boundary layers. 

4. NONPARALLEL STABILITY 

We confine our consideration to mean flows that are slightly nonparallel, that 
is, the normal velocity component V is small compared with the other components U m 
and W^. This in turn implies that all the mean-flow variables must be weak functions 
of the streamwise and spanwise coodinates. In other words 


U m ’ U m< x .- > - z ,> 




m 


m 


ITT 


(4.1) 


p m = 


P.n< x i* z i> 


T m 


yw’ z i) 


(4.2) 


where V* = 0(1). We describe the relatively slow variations of the mean-flow 
quantities in the streamwise and spanwise directions by the slow scales x x = ex and 
Zj = ez, respectively, where e = 1/R. 

We use the method of multiple scales (refs. 37,38) to determine a uniformly 
valid asymptotic expansion of the solution of equations (2.3)-(2.9), (2.12) and 

(2.13) in the form 


q(x,y,z,t;e) = [q Q (x 1 ,y,z 1 ,t x ) + eq^x^y^.tj) + — ]exp(i e) (4.3) 

where t 2 = et, a slow time scale, and 


30 

ax 


0 t(x ! ,Z X ) 


— = efx z ) — - - 


0) 


Assuming the phase e to be twice continuously differentiable, we have 


(4.4) 


3a _ 36 
3Zj 3X t 


(4.5) 
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In the case of mean flows that are independent of z, such as 20 mean flows and flows 
over infinite-span swept wings, a = a(x t ) and 6 = 6^), because the coefficients of 
equations (2. 3)-(2.9) are independent of Zj. Then, it follows from equation (4.5) 
that b = constant. Substituting equations (4. l)-(4.4) into equations (2.3)-(2.9), 
(2.12) and (2.13) and equating coefficients of like powers of e, we obtain problems 
describing the zeroth- and first-order disturbances. 

4.1 Zeroth-Order Problem 


p o T o, 


L , (u , v , w ,p ,T ) = - ip_(u) - all - sW ) (— - y - ) + ip_(au + Sw ) 
1 ' o o o o o' m v m trr v p m T m v o o' 


m m 


+ D(p v ) = 0 
v m o 


(4.6) 


L,(u ,v ,w ,p ,T ) = - ip m (w - all - eW )u + p V 0U m + i ap 
2V o o o o o' m' m m' o m o m o 


- 5 I’ + B X + 0(l, n, Du o> + W 1 + m > Dv o 

+ i av Du - a8y (1 + m)w + D(-r^ T ) DU + T D U ! = 0 
o m m v o x dT o' m dl o nr 


(4.7) 


.W,.p .T ) = - ip m (u> - all - sW ) v + Dp rt 
o o o r o o ' nr m nr o K o 

- | {ln m (l + m)(aDu 0 + bDw 0 ) + im(au 0 + BW Q )Du m 

- (“ 2 + ^)Vo + °(V Dv o) + 1 ( aDU m + sDW m) 3T V = 0 

(4.8) 


L,(u ,v ,w ,p ,T ) = - ip m (u - all - bW W + p v DW + isp r 
4 v o o o o o m v m m o m o m c 


- 1 (- + "> u 0 * ”)° v o + l8v o°"m 

- “ m (» 2 + r8 b"o + D <“ m D *'o> + “(ar T o)°V ar T o D \> * 0 


m o' 


(4.9) 
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co r 0 


MVVVP 0 - T o> = " V" - “ U m - sM m> T o + Vo 0T . + 

' <» - “ U » - BW J - R {^ m ( DU n. 0u o + “W 

♦ 2 V«° U * * sD W,)v 0 ♦ ^ [ (Du/ + (0« m ) ! ]T o } 

- f?l 2 V T o + u m [) T o - + 6 ^ ^ u m T o + 0 “in T ol ’ 0 


(4.10) 


u o = v o = w o = T o = 0 at y = 0 ( 4 - 11 ) 
u o ,v o ,w o ,T o - 0 as y - « (4.12) 
The solution of the zeroth-order problem can be expressed as 


u _ , Du. , v„ , T , DT , w„ , Dw rt 
L 0 0 0 0 0 0 0 J 


A(x 1 ,z 1 ,t 1 )c(x 1 ,z 1 ,y) 


(4.13) 


where c, is a column vector having 8 components and is governed by the quasi-parallel 
eigenvalue problem (3.13). The function A(x 1 ,z 1 ,t 1 ) is arbitrary at this level of 
approximation; it will be determined by imposing the solvability condition at the 
next level of approximation. 

4.2 First-Order Problem 


Using equation (4.13), we write the first-order problem as 
8 


Dz 


. - I a z. = D + E — + F — + GA 

In ^ nm lm n at L n ax x n az 1 n 

: 11 = Z 13 = Z 17 = 0 > Z 15 = 0 at y = 0 


3A 


3A 


(4.14) 

(4.15) 


z. is bounded as y -*■ « 
in 

where 


z. 


i 


' 1 2 


= Du 


i * 


' 1 3 


1 4 


= Pi 


(4.16) 

(4.17) 


' 1 5 


= T, 


' 1 6 


= DT, 


1 7 


= w, 


' 1 8 


= Dw, 


(4.18) 


and the D n , E n , F n , and G n are known functions of the c , a, 6, and the mean-flow 
quantities; they are defined in Appendix B. Since the homogeneous problem (4.14)- 
(4.16) has a nontrivial solution, the inhomogeneous problem has a solution only if 
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the right-hand side of equation (4.14) is orthogonal to every solution of the adjoint 
homogeneous problem (ref. 38) 


Dsi + l = 0 


m=l 


mn m 


r* = r* = r* = 0 r* 

**2 ^4 ’ 8 U » ‘*5 


= 0 at y = 0 


c* is bounded as y - <*> 
n 


The solvability condition takes the form 

J oJl I°nlf: +E nlf7 + F n ll7 + Wv 


(4.19) 

(4.20) 

(4.21) 


(4.22) 


Substituting for the D , E n , F^, and G n from Appendix B into equation (4.22) yields 
the following partial-cnfferential equation governing the modulation of A: 


aA . aA ^ aA 

— T~ + 00 -r CO 

at a ax b az 


ehj (a,8,X,z)A 


(4.23) 


where co^ and <o g are the group velocity components in the x and z directions; they, 
along with h lt are given in quadratures, as in Appendix C. 

Equation (4.23) shows that A is a function of x and z as well as a and 8. To 
determine the equations describing a and e, we differentiate equations (3.13) with 
respect to Xj and obtain 


3? n 8 ac 

ir) - l a 

3x i m=l 

3?! at 3 3c 5 3c 7 


8 aa. 


nf n i v , ( m-> _ ,-r- 3a , • c as r mn 

° 3x i " m i a nm^ax ^ lE n ax lF n ax, L g x , s 

1 m=l 1 1 1 m=l • a,8 


ax 

35 

ax, 


n . 


ax 


ax 


ax , 


= 0 at 


is bounded as y - - 


Imposing the solvability condition on equations (4.24)-(4.26) yields 


(4.24) 

(4.25) 

(4.26) 


CO 


9a 

ax, 


+ to . 


36 


= h. 


(4.27) 


where h 2 is defined in Appendix C. Similarly, differentiating equations (3.13) with 
respect to z y and imposing the solvability condition for the resulting problem, we 
obtain 


U) 

a 


3a_ a_6_ 

azj “e 3z t 


h 3 (a,6,X 1 ,Z 1 ) 


(4.28) 
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where h is defined in Appendix C. Using equation (4.5), we rewrite equations (4.27) 
and (4.28) as 


3a 3a . 

a 3X 8 3Z 2 



(4.29) 

38 , 36 _ h 

a 3X 8 3Z 3 



(4.30) 

Therefore, equations (4.23), (4.29), and (4.30) 

A, a, and 6 with x, z, and t. 

describe 

the 

modulation of 

For a monochromatic wave, 3A/3t = 0 and equation 

(4.23) reduces to 


3 A sA _ . a 

W a 17 + “6 Iz ~ £h > A 



(4.31) 

Nayfeh (ref. 21) argued that, in general, w /w is complex and hence equation (4.31) 
is elliptic for real x and z. In orcfer a that equation (4.31) be hyperbolic 

representing a propagating wave, wg/w must be real and equations (4.29)-(4.31) 
reduce to the ordinary differential equations 

dA _ ehlA da _ eh 2 de _ eh 3 

ds aj’dS ui * ds w 

a a a 



(4.32) 

along the characteristic 




dX 1 Anti dz 

ds " 1 and ds “ 

Thus, a 



(4.33) 

A = A 0 exp[ e J(h 1 /u) c( )dsj = A Q exp[ e J(h ,/^Jdx ] 
and 



(4.34) 

a = efh 2 ds, 8 = e/h 3 dS 



(4.35) 


where A Q is a constant. Therefore, to the first approximation, 

u h j 

u » A c (x,y,z)exp[ i f(a + 8 — )dx + ef — dx - iwt] (4.36) 

O 1 U) oj 

a a 

For the case of parallel flows, the condition wg/u) be real reduces to da/de being 
real, which was obtained by Nayfeh (ref. 21) and Tebeci and Stewartson (ref. 44) 
using the saddle-point method. 

4.3 Growth Rate 

Defining the growth rate in the streamwise direction as 


°x - Rea, l^ o™>] 
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we find from equation (4.36) that 


- - <«i + »i + e Rea ' <!r> + < 1n ^) 

a a 


(4.37) 


The first term is the quasiparallel contribution and the last two terms are due to 
the growth of the boundary layer. Whereas the second term is independent of the 
transverse direction, the third term strongly depends on the transverse direction. 
Hence, in contrast with the parallel or quasiparallel case, the local growth rate 
depends on both the streamwise and transverse directions. Moreover, the transverse 
variation depends on the flow variable being investigated. 


It follows from equation (4.36) that the amplification factor is 


c n (x,y,z) r x o> x h L 

, ■ fx y z) exp l ’ -f (“i + 6 i ^) dx + eReal J ~ dx l ( 4 - 38 ) 

W y,Z ' X 1 "’a X U a 

0 0 

Again, the amplification factor is a function of the transverse location. However, 
the variation is a simple one. For a detailed discussion of the case of a wave 
packet, the reader is referred to the article of Nayfeh (ref. 21). 


For the case of two-dimensional flows, 8 is a constant. Then, for the case of 
monochromatic waves, equation (4.36) reduces to 


u « A o c 1 (x,y)exp[i(Jo(dx + ez 


U)t ) 



(4.39) 


Consequently, the growth direction is the streamwise direction and the growth rate 


a = Real [f^(lnu) ] 

is given by 


O = - a. + eReal(^) + Real [|^( Ins : ) ] (4.40) 
a 

Since is a function of y and, in general, distorts with the streamwise 
distance, one may term stable disturbances unstable and vice versa. Because of the 
mode-shape distortion, neutral stability points are a function of both the transverse 
and streamwise positions. The experiments of Laufer and Vrebalovich (ref. 23) 
clearly show that the growth of the disturbances and the neutral curves obtained 
depend on the value of y*/L* at which the observations were made. Moreover, a 
different growth rate would be obtained if one replaces u with another variable such 
as v, p, w, or T. On the other hand, for the case of parallel flows, the last two 
terms in equation (4.40) vanish and the growth rate is unique and independent of the 
varible being used or the transverse direction at which the growth rate is 
determined. Therefore, to compare the analytical results with experimental data in 
growing boundary layers, one needs to make the calculations in the same manner in 
which the measurements are taken. Available experimental stability studies almost 
exclusively use hot-wire or hot-film anemometers following disturbances in the 
boundary layer. The hot-wire or hot-film response is a combination of velocity. 
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density, and temperature fluctuations. The hot-wire response when operated at high 
constant overheat is proportional to the rms of the mass-flow fluctuations. To 
describe fully the disturbance, one needs measurements of all fluctuation 
characteristics such as the rms amplitude, spectra, wave angle, and propagation speed 
as functions of x, y, and z. Except for the experiment of Kendall (ref. 25), no 
information is available on the disturbance waveangle. Laufer and Vrebalovich 
reported measurements at different constant y*/L*, whereas Kendall, Demetriades, 
Legiga et al, and Stetson et al reported measurements at one constant y*/L* located 
in the wideband energy peak. However, available calculations do not agree with each 
other. Some calculations show large nonparallism influence, whereas other 

calculations show small influence. 

5. SUBHARMONIC INSTABILITY 

In this case, the basic flow is taken as the sum of the mean steady flow and a 
two-dimensional quasiparallel T-S wave; that is 


u b 

■ V*’ 

+ A[ tl (y)e 10 + cc] 

(5.1) 

v b 

= A[ c 3 (y)e ie + cc] 

(5.2) 

p b 

= p m + 

A[t 4 (y)e ie + cc] 

(5.3) 

T b 

- v*> 

+ A U 5 (y)e 16 + cc] 

(5.4) 

w b 

= 0 


(5.5) 

30 
3 X 

°v 

and -|f- = - a) 

3t 

(5.6) 


where cc stands for the complex conjugate of the preceding terms, a and u> are real, 
and A and a are approximated locally by constant values. Substituting equations 
(5. l)-(5.6) into equations (2.3)-(2.9), (2.12) and (2.13) yields a system of partial 
differential equations whose coefficients are independent of z, periodic in x and t, 
and dependent in a complicated manner on y. Consequently, the z-variation can be 
separated and using Floquet theory, one can represent the solutions of the problem as 


u = exp(o x x + o t t)cosez ^(x.y.t) (5.7) 
v = exp(o x x + a^t)cosez i> 3 (x,y,t) (5.8) 
p = exp(o x x + o t t)cosBZ <Mx,y,z) (5.9) 
T = exp(o x x + o^t ) cosbz i> 5 (x,y,t) (5.10) 
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(5.11) 


w = exp(a x x + a t t)sinez <t> 7 (x,y,t) 

where a and a. are called the characteristic exponents and the <t> are periodic in x 
and t. x For tne subharmonic parametric case, the <j> n have a perioa that is twice the 
period of the primary flow and to the first approximation equations (5.7)-(5.11) 
become 


exp(a x x + a t t)cosBz[c 1 (y)exp(^ie) + cc] 

(5.12) 

exp(a x x + o t t)cos6z[c 3 (y)exp(-^i0) + cc] 

(5.13) 

exp(a x x + a t t)cos8z[s 4 (y)exp (-^ie) + cc] 

(5.14) 

exp(a x x + 0 t t)cosez[s 5 (y)exp(-^ie) + cc] 

(5.15) 

exp(o x x + a t t)sinez[s 7 (y)exp(-^ie) + cc] 

(5.16) 

case of temporal stability a = 0 and o. * 0, 

.. _ n i n X L 

whereas for the case of spatial 


stability = 0 and o x * 0. 


Substituting equations (5.1)-(5.6) and (5. 12)-(5. 16), into equations (2.3)-(2.9), 
(2.12) and (2.13) and equating the coefficients of exp(^-ie) on both sides leads to 
the following problem governing the s and either a or o C . (ref. 45): 

is A U 


l°t ~ l lu) + ^°x + l ia 

Dp 

+ — C 3 + Dc 3 + 8^ 7 
M m 

+ (° x + | ia r ) A ^ i(y m J; 

+ - Vs)l 


r> U ,J( YM X " p m^s) + K + I ia r )f ^ 

+ (o x + - P^s)^! 

** " p m^) + p^ ay l p n/ YM ® C4 " p m ? s)^3 

+ a (^ 4 - P m 4 5 )€ 7 = 0 


(5.17) 
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I I I I ^ 

K ■ + U m (a x + I 1o r) - 7% (a x + i ia r )2 + fr ki 

m m 


, p m . . r pm y m f , 1 . ,| r , r 

lDU m ' ^TR (o x + 2 1a r } ^3 + 7 ~ ^ 


D(y;OU) 

'mm' 


m 


m 


p m R 

m 


u'DIL (m+l)u m s 
m m v ' m 

pR m . 

m m 


1 1 
p R < a x + 2 ia r^ “ p R K + 2 ia r )D5 3 


m 


JJ 4 4 

- ~Tr D? 2 + A ^i(°x + 2 1a r^ + K + 2 i(i) + U m (a x 

m 

- |i“ r )l (^ 4 - P m ^ 5 ) - (a* + 

m 

D ( y'C 5 ) _ 2 

* # ['3 - -fr ^ + + DU m (,M. t> - „ m c s ) 

m 

d (^ 5 ) , i. 


m 


R K ~ 2 iot r^ 5 3 + A t i ^°‘r U m " “) c i 


2 

+ d V 3 H yM J 4 - P m ^ 5 ) - [Pm (ir Vi + mD ^)(° x + I ia r ) 

Ay' 
tr 

j_F 
m 


A 


1 . 


m 


+ p m^ ^ 1 + l 0 l r^ c 3 ) + + ^“r^ 3 ) Us p r (^Ci 


Agy. 


+ ia^He ~ 7~R 21 [ m ( a x + i ia r ) + (° X ‘ ^<V)k 7 = 0 


1 , 


i ** 3 ' 4 6 p^R 


I 

i 


I 


(5.18) 
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2( y -1)M 2 DU 2( y -1)M 2 DU . 

® m _ , r nT v ' co m, . li \i r 

^ DT tn o_R x 2 r^ c 3 
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(y-1)M 

p m 


pR 
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K ‘ 1 it0 + U m (o x + 1 iot r)^ + W ' 1 


1 0) 


+ U m (o x + 2 ia r ) " 
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m m mm m J J 5 
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0 RPr p RPr ^6 
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05, 


(y- 1)M 2mv , 

+ A t i Vs — [ ia r ^4 + nr 1 (° x ■ 2 io v^ 3 
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2iry a„ , 2 A(y-1)M 

R (°x " 2 " Rp t Dc > + ia r c 


m 


r w 3 


2 2 ( Y - 1 ) M , 

+ ^m DU m^5)52 + A ((^ 4 - P m ^)DT m + Dc 5 - (° x - 2 


m 


+ iot r C 3 + ^s DU JH 3 + A{ Y [i(ot r U m - u)c 5 3- C 3 DT m ] 


(y-1) 


" — (° x - \ ia r )^i^4 + M* yM^[ i (ot U - u»)c 5 + c 3 DT n 

p m 


1 1 , 2( y - 1)M p'DU 

+ [o t + 2 1w + u m (°x - 2 la r ) !(y m ^4 - P m ^s) - — 


Rp 


(Oc L 


Dp p 2 1 2 1 

+ ia r^) - Dc s - U° x + 4 a r } ^ + ia r (a x + 2 ia r )c 5 

m m 

+ o\ 5 - 6 2 ; 5 ]}^ 5 + A{^3 - -jp- [(D(p'^ 5 ) + u'Dc 5 )]}c 6 


2 m 2 

2m( Y -l)M Su m 2 ( y - 1)M u 

A p R ( ia r £i + °^ 3 )^7 " A T“r ( rD ^3 + 1nta r c l )D5 3 
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A(y-1)M‘ 








°m " p m RPr 
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= 5 3 

= C 7 = 0 and ? 5 

= 0 at y = 0 

(5.22) 

^n 

* 0 

as y ® 


(5.23) 

where 

^ 2 

= De, 

l, 5 6 = D^ 5 , and 

’ 8 = D ^7 

(5.24) 

As 

an 

alternative to the 

collocation technique 

used by Herbert for the 


incompressible case (refs. 46,47), Nayfeh and Harper (ref. 45) followed Nayfeh and 
Masad (ref. 48) and used a shooting technique to solve the eigenvalue problem 
consisting of equations (5.17)-(5.23) after casting them into a first-order system of 
equations. The boundary conditions (5.23) at infinity were replaced with boundary 

conditions at some finite value y ma of y outside the boundary layer. Thus, for y > 

y max , U -*• 1 and DU - 0. Since tne primary eigenfunctions c decay exponentially 
outside m the boundary layer, Nayfeh and Harper (ref. 45) followed Nayfeh and Masad 
(ref. 48) and chose a large enough y m so that the t are very small. The error 
involved was monitored by choosing a certain value of y™ x , and solving the problem, 
then increasing y_ ax , resolving the problem and noting tne effect of increasing y max 
on the accuracy of the solution. 

Since no numerical or experimental results are available for compressible flows, 
Nayfeh and Harper checked the compressible code against the incompressible results of 
Herbert (refs. 46,47) obtained using a collocation method and those of Nayfeh and 
Masad (ref. 48) obtained using a shooting technique, which are in good agreement with 
the experimental results of Kachanov and Levchenko (ref. 49). The agreement is 

excellent. Next, Nayfeh and Harper produced results to evaluate the influence of 

compressibility on the secondary instability of first- and second-mode waves. 

Letting R = 1048.8 and choosing an F = 83xl0 -6 , Nayfeh and Harper solved the 
eigenvalue problem and obtained the 2D T-S waves for M^ = 0.01, 0.8, and 1.6. These 
primary waves are stable having the spatial decay°°rates -0.0216, _-0.0166, and 
-0.0066. Setting the amplitudes of these primary waves at a = 0.02//2, choosing a 
value for 6, and assuming a value for y., they solved the secondary eigenvalue 
problem and obtained y . Varying y. and repeating the calculations, they determined 
the maximum value of y r over all possible values of y.. We note that when M > 0 
the maximum value of*" y does not correspond to yl = 0 as in the incompressible 
case. The results are Shown as the dark curves ir\ Figure 15. The growth rates 
corresponding to a = 0.01//2 are shown as the light curves in Figure 15. Figure 16 
shows the growth rates when a = 0, that is, the growth rates of the free propagating 

waves. The results in Figures 15 and 16 show that the secondary growth rates are 

much larger than the growth rates of the free waves. In fact, for M = 1.6, the free 
wave is stable for all values of e, whereas the secondary wave is unstable for a very 
wide range of values of s. Increasing the amplitude of the primary wave results in 
an increase in the growth rates of all secondary waves. For a given Mach number, the 

growth rate has a broad maximum. This maximum shifts to larger values of 6 as the 

amplitude of the primary wave increases. And for a given amplitude of the primary 
wave, the maximum growth rate shifts to lower values of 6 as the Mach number 
increases. 

As discussed in Section 3, as M increases beyond 3.0, higher-modes become 
unstable and dominate the instability at large Mach numbers. Figure 17 shows the 
variation of the growth rates of secondary waves with spanwise wavenumber when the 
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primary wave is a first- or second-mode wave. In both cases, F = 120xl0~ 6 , a = 
0.01//2, and y. = 0. The first-mode primary wave and its corresponding secondary 
waves are calculated at R = 1150; the primary wave is slightly unstable. The second- 
mode primary wave and its corresponding secondary waves are calculated at R = 1950; 
the primary wave is very stable. It follows from Figure 17 that the maxima of the 
growth rates of the secondary waves are comparable. Figure 18 shows the growth rates 
of the free waves (secondary waves when a = 0). The free wave corresponding to the 
second-mode wave is stable for a wide range of s, whereas the free wave corresponding 
to the first-mode wave is slightly unstable for the same range of 8. Comparing 
Figure^ 17 and 18, we conclude that a primary first-mode wave having an amplitude of 
0.01//2 increases the growth rate of the secondary wave by an order of magnitude. 
Moreover, the presence of a primary second-mode wave strongly destabilizes the 
secondary wave. 

When the primary wave is a first-mode wave it follows from Figure 17 that the 
maximum growth rate of the secondary wave occurs when g = 0.199. Figure 19 shows the 
variation of the secondary growth rate with frequency. The maximum growth rate 
increases by about 25 % and shifts to an F « 112x10 . 

6. EFFECT OF IMPERFECTIONS ON STABILITY OF FLOWS OVER PLATES 

The boundary layers over natural laminar-flow components in the presence of 
surface imperfections (e.g., waviness and steps) must accurately be computed so that 
the effect of these imperfections on the stability and transition can be evaluated. 
Moreover, the magnitudes of the imperfections under consideration are such that 
strong viscid-inviscid interaction and small separation bubbles are unavoidable. 
Definitely, solutions to the full Navier-Stokes equations can accurately predict such 
flowfields provided that the grid is so fine that important flow structures are not 
smeared by the truncation errors or artificial dissipation. However, the number of 
flow cases that needs to be investigated is very large, and this makes solving the 
full Navier-Stokes equations a very expensive task. A more economical alternative is 
to solve the interacting boundary- layer equations. 

Ragab, Nayfeh and Krishna (ref. 50) investigated the accuracy of the 

compressible interacting boundary- layer computations and their limitations in 
predicting flows over surface imperfections. They compared the results of 
interacting boundary-layer computations with solutions to the full Navier-Stokes 
equations. Comparisons were made for the mean flow profiles as well as the stability 
characteristics such as the growth rates and amplification factors of linear 
stability waves. 

The thin-layer compressible Navier-Stokes equations are solved using the well 
known computer code "ARC2D" which has been developed at NASA Ames (version 1.5 GAMMA, 
72/85). The code incorporates different methods of solutions all of which are 
implicit in time, and it uses second-order central differences in space. Ragab et al 
(ref. 50) selected the method of solution in which the diagonal form of the equations 
are used (refs. 51,52). Mixed second- and fourth-order dissipation terms were added 
explicitly and implicitly, and the obtained pentadiagonal system of equations was 
solved directly. 

Sheared Cartesian grids were used for all the cases presented here. An example 
is shown in Figure 20 for a smooth backward-facing step. The equation of the step is 

y = j h(l + erf x), x = Re _3 ^ 8 x 5 ^ 4 (x - 1) (6.1) 
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where Re is the Reynolds number based on the distance from the leading edge to the 
step center (x = 1) and x = 0.332057. The step specified by equation (6.1) is 
originally given by Smith and Merkin (ref. 53) who analyzed the incompressible flow 
using triple-deck theory. The numerical values used in Figure 20 are Re = 10 and h 

= - 0.003. We note that the y-coordinates in Figure 20 are magnified by a factor of 

20 relative to the x-coordinates. The inflow boundary of the computational domain is 
at x « - 0.06 (i.e., the plate leading edge x = 0 is included in the domain), and the 
outflow boundary is located at x = 2.0. The top boundary is placed at y = 0.4. More 

details about the grid will be given in the results section. 

Ragab, Nayfeh and Krishna developed a code for solving the compressible 
interacting boundary- layer equations. The numerical method is an extension of 
Veldman's method (ref. 54) to compressible flows. The salient feature of the method 
is the simultaneous solution of the boundary- layer equations and the inviscid flow, 
which is given by the small disturbance theory of compressible potential flow. A 
similar treatment is presented by Davis (ref. 55) and Nayfeh, Ragab, and Al-Maaitah 
(ref. 56). 

In this subsection we compare the mean flow predicted by the Navier-Stokes code 
ARC2D and the interacting boundary- layer code IBL for the flow over a backward-facing 
step. The step is specified by equation (6.1). The 6 step height is h = - 0.003, the 
Mach number is 0.5, and the Reynolds number is 10 . The wall is assumed to be 
insulated. 

The grid used in the IBL calculations has a uniform streamwise step size of ax = 
0.005 and a geometrically stretched grid in the n- direction, where n is the Levy- 
Lees variable. At the wall An = 0.05 and the stretching factor is 1.05. Four grids 
are used with the Navier-Stokes solver. In grid 1 (136x70), a uniform ax = 0.005 is 
used in the range 0.9 < x < 1.1. For x < 0.9 and x > 1.1, ax is stretched 
geometrically at the rate of 1.05 provided that ax does not exceed 0.03. If ax 
exceeds 0.03, a uniform spacing ax = 0.03 is used. The step size in the y direction 
is geometrically stretched between the wall and the top boundary y « 0.4 with Ay. = 
1.5x10" . In grid 2 (136x99), the streamwise grid is the same as in gr^d 1, and the y 

grid has 99 points between the wall and y » 0.4 with Ayj = 0.3x10“ . In grid 3 

(166x120), the streamwise grid in the interval 0.9 < x < 2.0 is the same as that in 
grid 1 and more points are added in -.06 < x < 0.9 so that the number of streamwise 
points is 166, and the y grid has 120 points between the wall and y « 0.4 with Ay^ = 

0.3xl0"\ In grid 4 (176x153), the streamwise grid is the same as that in grid 1 

except ax = 0,003, and the y grid has 153 points between the wall and y « 0.8 with 
= 0.3x10 . 

The friction-coefficient distribution on the step surface is depicted in Figure 
21. We observe that there is an appreciable difference between the results predicted 
by ARC2D using the different grids. The trends in Figure 21 suggest that the 
agreement can be improved slightly if a finer grid is used with ARC2D. The results 
obtained using the IBL code are in good agreement with those obtained using ARC2D 
with the finest grid. Figure 22 compares the pressure-coefficient distribution 
predicted using the IBL code with that obtained using ARC20 with the finest grid. 
The agreement is fairly good. Thus, we conclude that the IBL code can be used to 
calculate accurately the mean flows over steps, humps, and troughs as long as massive 
separation and vortex shedding do not occur. 

The stability calculations demand highly accurate velocity and temperature 
profiles - their magnitudes as well as their derivatives. Figure 23 compares the 
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spatial growth rates -a. calculated at x = 0.64 for an F = 80x10" using the I BL 
predicted profiles with chose calculated using the ARC2D predicted profiles with the 
finest grid. The agreement is very good, indicating that the IBL code predicts 
fairly well the magnitudes as well as the derivatives of the velocity and temperature 
profiles. Using these growth rates, we calculate the amplification factor N from 

R 

N = - f 2a. dR 

D 1 


where R = /Re . The variation of the N factor calculated using the growth rates 
based on the profiles obtained with the IBL code is compared with those calculated 
using the growth rates based on the profiles obtained with the ARC2D code with the 
four grids in Figure 24. The latter N factor distribution seems to converge as the 
grid is refined to a distribution that comes closer and closer to that obtained using 
the growth rates based on the profiles predicted by the IBL code, providing a further 
credibility to the IBL procedure. 

Figure 25 shows the influence of the Mach number on the N-factor distribution 
for the most amplified 2D wave for a backward-facing step having an h = - 0.003 and a 
slope of -4.349°. It is clear that compressibility is stabilizing. 

In Figure 26, for = 0.5, we show the influence of the disturbance frequency 
on the N-factor distribution for waves propagating past the same backward-facing 
step. It is clear that the most dangerous frequency F is about 50x10" . 

7. SWEPT WING BOUNDARY LAYERS 

All of the numerical examples discussed in the preceding sections are concerned 
with the 2D and 3D stability of 2D or ax i symmetric boundary- layer flows. In this 
section, we discuss the stability of 3D boundary layers and in particular boundary 
layers on transonic sweptback wings because of their aeronautical importance. The 
stability of 3D boundary layers differs from that of 2D boundary layers in that a 3D 
boundary layer is subject to crossflow instability and in that the growth direction 
in a 3D boundary layer need not coincide with the wave direction. 

From his flight tests on aircraft with sweptback wings, Gray discovered that the 
boundary layer became turbulent closer to the leading edge than on a corresponding 
unswept wing. Using evaporation methods for indicating the state of the boundary 
layer, he discovered the existence of regularly spaced vortices whose axes lie in the 
streamwise direction. Since on a swept wing the spanwise pressure gradient deflects 
the boundary layer toward the region of low static pressure, the flow paths of the 
boundary- layer profiles differ from the potential flow streamlines and a crossflow 
developes in the direction normal to the streamlines (i.e., the mean flow is three- 
dimensional). The crossflow profiles have inflection points, making them dynamically 
unstable and hence leading to the generation of the vortices. This crossflow 
instability was confirmed in wind tunnels on large swept wings by Gregory, Stuart and 
Walker (ref. 34). 

The feasibility of using suction to maintain laminar flow in the presence of 
crossflow instability was shown by Pfenninger et al (ref. 57), Bacon et al (ref. 58), 
Gault, and Pfenninger and Bacon (ref. 59) on a 30-deg sweptback wing. This 
feasibility culminated in the successful maintenance of full-chord laminar flow on an 
X-21 wing. 
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The problem of laminar-flow control on sweptback wings and the discovery of Gray 
of the crossflow instability stimulated research into the linear stability of three- 
dimensional flows. In contrast with the problem of 2-D flows, disturbances in a 3-D 
flow are always three-dimensional. Stuart (ref. 34) derived the general linear 
equations that describe the 3-D stability of 3-D incompressible boundary layers over 
bodies, including the effects of boundary-layer growth and body and streamline 
curvatures. He showed that the temporal parallel stability problem of a 3-D 
incompressible flow can be reduced to that of a 2-D stability problem. For a given R 
and u>, his transformation relates the wave numbers a and e in a 3-D flow with the 
strgamwise and spanwise velocity components U m and W ffl to the wavenumber k = (a 2 + 
s 2 ) ^ in a 2-D flow having the velocity U = (aU m + BW m )/k. He used the inviscid form 
of these disturbance equations to explain the crossflow instability as an 
inflectional instability. Brown (refs. 7,16) numerically solved the temporal viscous 
eigenvalue problems for flows over a rotating disk and a few sweptback wing boundary 
layers with distributed suction. 

Since the uncoupling of the crosswise stability problem from the streamwise 
stability problem is artificial, a number of attempts have been directed toward 
determining the stability of the combined streamwise and spanwise stability problems 
to determine the suction requirements for maintaining laminar flow over sweptback 
wings. The form of the 3-D disturbance is given by equations (4.3) and (4.4) leading 
to the dispersion relation u> = w(ct,8,x,z) . We note that, in contrast with the case 
of parallel flows, a and 6 are functions of x and z and the phase function e cannot 
be written as ax + sz - tot. Instead a, e, and u are related to e as in equations 
(4.4). 

To complete the problem formulation, one needs to specify initial conditions in 
addition to the governing equations and boundary conditions. For the case of a 
general pressure disturbance at say the curve x = a. 


P(x = a, y,z,t) = P(y,z,t) (7.1) 

This arbitrary disturbance will generate a continuum of wave components of the form 
given by equations (4.3) and (4.4) and the observed wave motion consists of their 
superpositions. For the case of spatial stability in parallel mean flows, the 
pressure disturbance can be expressesd 


p(x,y,z,t) = JJ c 4 (y;8,w)exp[ia(8,w)x + iez - i<i)t]du>de 


(7.2) 


On the other hand, if the initial condiitons are generated by a source (such as a 
vibrating ribbon) oscillating at the frequency w (i.e., monochromatic wave) at the 
curve defined by x = a, then 


P (x = a,y,z,t) = P(y,z)exp(-iu>t) (7.3) 

Then, for the case of parallel mean flows, the pressure disturbance can be expressed 
as 


p(x,y,z,t) = exp(-icot)f Jc 4 (y;e)exp[ ia(e)x + iezjdB 


(7.4) 
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To evaluate the integral in equation (7.2), one needs to solve the eigenvalue problem 
for all possible values of 6 and u> and then perform the integration, an expensive 
procedure. Similarly, to evaluate the integral in equation (7.4), one needs to solve 
the eigenvalue problem for all possible values of e and then perform the integration, 
again an expensive procedure. However, if the disturbances are sufficiently weak 
they will not influence transition until they have traveled a large distance from the 
source. Then, asymptotic methods, namely the saddle-point method, can be effectively 
used to determine the following approximations to the integrals in equations (7.2)and 
(7.4): 


2 A 

l ae3a) J J 


3 a 3 a 
X l 2 2 


where 


z 

X 


30 3 a) 

3a , t 3a 

36 3nd X = 1Z 


and 


P(x,y,z,t). explM^x.^z-u-t 


3 2 o t 

as' 


( 6 *) 


where 


z _ d<* 

x ds 


(7.5) 


(7.6) 

(7.7) 


(7.8) 


It is clear from equations (7.6) and (7.8) that for a physical wave, 3a/38 and 3a/aw 
must be real for an arbitrary disturbance and da/ds must be real for a monochromatic 
disturbance. 


The preceding discussion shows that by the time the wave motion is important for 
transition, it has already evolved to either the state given by equation (7.5) in the 
case of an arbitrary disturbance (corresponding to natural transition) or the state 
given by equation (7.7) in the case of a monochromatic disturbance. In both cases, 
the disturbance is dominated by a single wave component having a fixed w, e, and a. 
Using the method of multiple scales or the method of averaging to account for a weak 
growth of the boundary layer in x and z, one finds that by the time the wave motion 
is important for transition (i.e., for large x), it has evolved again to the state 
given by equations (4.3) and (4.4); that is, the disturbance is dominated by a wave 
packet centered at the frequency w and spanwise wavenumber b. However, in this case, 
a, 6, and A are not constants; equations (4.23), (4.29) and (4.30) describe their 
modulation with x, z, and t. To solve these modulation equations (Cauchy problem), 
one needs to specify initial conditions on non-characteristic curves. 

If the initial data consist of a wavepacket centered at the frequency w, then 


A (x = a,z,t) = A 0 (z,t), 8 (x = a,z) = s 0 (z) (7.9) 

and it follows from equations (4.23), (4.29) and (4.30) that the wave propagates 
along the characteristics 


dt i dx dz 

ds " 1# ds V ds w e 


(7.10) 
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For a physical problem. 


a) and (a) must be real. 

a 8 


Along these characteristics. 


dA _ . . 

ds ‘ eh i A ’ 


da 

ds 


eh, , 



(7.11) 


If the initial data consist of a single frequency (i.e., monochromatic wave), then 


A(x = a ,z) = A Q (z), s(x = a,z) = s Q (z) (7.12) 

It follows from equations (4.29)-(4.31) that the wave propagates along the 
characteristics defined by equations (4.33), where u) g /u) a is real. The modulation of 
A, a, and 8 along these characteristics is governed a by equations (4.32). The 
solution of equations (4.33) is a one-parameter family of curves. 

The question arises how one can use the results of the linear stability theory 
for transition prediction or evaluation of proposed laminar-flow concepts and 
requirement on swept wings? Motivated by the success of the so-called e N method in 
correlating experimental transition data in incompressible 20 flows, many 
investigators have attempted to extend it to 30 incompressible and compressible 
flows. Whereas in 2D flows, the direction of growth is known, namely, the mean flow 
direction, in 3D flows, the direction of growth is not known a priori. In 20 flows, 
one can either use spatial stability theory and calculate the spatial growth rate 
directly or use temporal stability theory, calculate the temporal growth rate, and 
then calculate the spatial growth rate using the group velocity. Integrating the 
spatial growth rates yields the N factors. 

Srokowski and Orszag (ref. 60) developed a computer code called SALLY using the 
temporal incompressible stability theory. At any location on the wing, the code 
iterates on a and 8 (which are taken to be real) to yield the maximum temporal growth 
rate for a given dimensional frequency m . Then, several frequencies are examined 
until a global local maximum is obtained. 1 " The resulting amplification rate is then 
converted into a spatial amplification rate using the real part of the group 
velocity, which is then integrated along a trajectory defined by the direction of the 
real part of the group velocity. Hefner and Bushnell (ref. 61) calibrated the 
results of the envelope method using available transition data on swept-back wings. 
Later, Malik (ref. 36) developed a code called COSAL, which is the compressible 
version of SALLY. Although the envelope method is easy to implement, it is 
artificial because one continuously hops from one wave to another. Moreover, the N 
factors calculated using the envelope method increase monotonically to the end of the 
instability region, and hence the envelope method may over estimate the amplitude 
ratios. 

For an infinite swept wing, the mean profiles are independent of the spanwise 
coordinate, and hence it follows from equations (4.4) and (4.5) that 8 is constant. 
Mack (ref. 62) performed crossflow stability calculations over the forward region of 
a 35° sweptback wing using this condition and determined the N factors for a zero 
frquency and a band of initial wave numbers. He found that the peak value maximized 
over all wave numbers is 7.8 compared with 11.2, obtained using the SALLY code. In 
contrast with the results of the envelope method, which increase monotonically in the 
instability region, the results of Mack achieve a peak because he followed a given 
wave from a more unstable region to a less unstable region. Using the temporal 
theory, Lekoudis (ref. 63) and Mack (ref. 64) determined the influence of suction on 
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the stability of flows over sweptback wings. Lekoudis (ref. 65) ascertained the 
influence of cooling on the stability of flows over sweptback wings. 

If the wing is not infinite, the mean flow will be a function of both spanwise 
and streamwise coordinates and hence a and 6 are functions of x and z. In this case, 
the wave motion will be dominated by a single frequency spatial wavepacket as 
discussed above. Based on this theory, Nayfeh (ref. 21) proposed a method for 
determining the most amplified disturbance propagating from some given initial chord 
location on a wing. We select one specific wave at an initial point and then follow 
only that one wave along its trajectory. Then, we change the initial conditions and 
repeat the calculations to determine the most unstable wave. 

To demonstrate the method we choose a*, the dimensional form of a r , to be a 
constant at some initial chordwise location x = a for all z. At this initial point we 
still have three unspecified values a^, a r , and c^. Two of these are determined 
through numerical integration of the disturbance equations and satisfaction of the 
boundary conditions through a Newton-Raphson iteration. The third is found by 
requiring that w A> be real. This fixes the direction of marching. From here we 
evaluate 6 a 


ds u> 

a 

ds = 

ds e a) 

a 

To stay on the initial wave, we increment the characteristic variable s by ds, 
evaluate dx and dz from equations (4.33), and correct a and 6 by 

^2 

a = a + e — ds 
0) 

a 
^3 

8 = 6 + e — ds 

U) 

a 

We then renondimensional ize with respect to local edge variables at the new x and z. 
To test that we are on the same wave, we integrate the disturbance equations at the 
new x and z with the new a and s to see if the boundary conditions are satisfied and 
u)g/u) a is real. If they are, we evaluate da/ds and ds/ds, increment s, correct a and 
6, and continue marching along the trajectory. If these conditions are not all 
satisfied, we decrease the step size ds until they are and then proceed as above. 

disturbance initially, it follows from equations 


(7.17) 

% 

s — ) + e — idx - iwt} (7.18) 

U) ' td ‘ 

a a 


Since A is the amplitude of th 
(4.32a) that 


Hence 


A = A 0 exp( J e — ds) 
at a 


u = AjjC^Xj.y.ZjJexp { J [i(a + 


and the n-f actor is given by 


(7.15) 

(7.16) 


(7.13) 

(7.14) 
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As we march we compute the n-factor from where the disturbance first goes unstable 
< t0 wh ere it becomes stable again < 0). Once we compute n we change e r 

i *r 

to determine which e is associated with the largest value of n. This tells us the 
initial spanwise wavenumber of the most dangerous wave for a given frequency from a 
given chord location. Then we vary the dimensional frequency to determine the most 
dangerous frequency. 


This approach has been implemented by Padhye and Nayfeh (ref. 66) and Reed and 
Nayfeh (ref. 67) for the X-21 wing with upper surface chordwise pressure coefficient 
and suction distribution as shown in Figures 1 and 2. A schematic of the wing is 
presented in Figure 3 showing the coordinate system used. This wing was designed for 
laminar-flow control and derived from the NACA 65A210 airfoil. The sweeps at the 
leading and trailing edges are 33.2 and 19.1 degrees, respectively, and the 
freestream velocity is 774.4 feet per second. Under the assumption of negligible 
spanwise (along roots) pressure gradient, there are no twist and tip and wing-body 
effects included in the analysis. Also they made the assumption of constant Prandtl 
number and specific heat at constant pressure. 

Reed and Nayfeh (ref. 67) focused their attention on the aft region of the wing, 
specifically the 70-percent location on the 14.66-foot chord. They introduced 
disturbances at this point and followed them along their trajectories as described 
above to ascertain their stability characteristics. This initial location lies in 
the rear adverse chordwise pressure gradient region. 

For a specific constant dimensional frequency introduced at the 70-percent 
chord, to find the spanwise wavenumber maximizing the amplification factor, they 
consider the wavenumber resulting in local maximum temporal growth rate. That is, 
they took the frequency w to be complex and the two wave numbers a and e to be 
real. They specified ay and iterated on a, e and to satisfy the boundary 
conditions and the condition that cy/w be real. Then they converted to spatial 
stability and marched along the direction defined by the group velocity, which is 
kept real as described above. The amplification factor was computed in marching. 
This process was done for a whole gamut of frequencies until the frequency with 
maximum amplification factor was identified. Figure 4 shows results of frequency 
plotted versus spanwise wavenumber causing maximum growth for disturbances introduced 
at the 70-percent-chord location. The distribution appears to be linear. 

As a result of searching, they found that the most unstable disturbance 
introduced at 70-percent chord is the one with a frequency of about 330 hertz and 
spanwise wavenumber of 272 per foot. Figure 5 shows disturbance amplification ratios 
for frequencies in this neighborhood. If the stability was performed using 
incompressible theory, the most unstable frequency would be about 200 Hz. The 
addition of compressibility increases the frequency of the most unstable 
disturbance. Moreover, compressibility significantly reduces the growth rates and 
amplification factors. 
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Figure 1. Multiple wave numbers of 2D noninflectional neutral waves (c=l). 

Insulated wall, wind tunnel temperatures, calculations of Mack (ref. 
11 ). 





Figure 2. Pressure-fluctuation eigenfunctions of first six modes of 2D 
noninflectional neutral waves (c=I) at M = 10. Insulated wall, T* 

= 50°K, calculations of Mack (ref. 11). 
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Figure 3. Multiple wave numbers of 2D inflectional neutral waves (c=c_). 

Insulated wall, wind tunnel temperatures, calculations of Mack (ref. 
11 ). 



Figure 4. Pressure-fluctuation eigenfunctions of first six modes of 2D 

inflectional neutral waves (c=c.) at M^ = 10. Insulated wall, T* = 

50°K, calculations of Mack (ref. II). 
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Figure 5. 


Effect of Mach number on maximum temporal amplification rate of 2D waves 
for first four modes. Insulated wall, wind tunnel temperatures (ref. 
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Figure 6. 


Temporal amplification rate of first and second modes versus frequency 
for several wave angles at = 4.5. Insulated wall, T* = 311°K (ref. 
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Figure 7. Effect of wall cooling on ratio of maximum temporal amplification rate 
with respect to both frequency and wave angle of first and second modes 
at M = 3.0, 4.5 and 4.8 to insulated-wall maximum amplification rate. 
Wind°°tunne1 temperatures (ref. 11). 



Figure 8. Two-dimensional neutral-stability curves at = 1.6, 2.2, 2.6, 3.0, and 
3.8; is the first neutral inviscid wavenumber with c = c $ (ref. 11). 
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Figure 9. Distribution of maximum temporal amplification rate with Reynolds number 
at (a) M m = 1.3, (b) M = 1.6, (c) M = 2.2 and (d) = 3.0 for 2D and 

3D waves? Insulated wall, wind tunne”) temperatures (ref. 15). 
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Figure 10. Variation of the maximum spatial amplification rate (with respect to 
frequency) with streamwise position at = 1.6 for an oblique wave 

at ii> = 50° and a two-dimensional wave (ref.” 19) 



.7 



Figure 11. Effect of wall cooling on the most amplified first and second mode 
disturbances in a flat plate boundary layer at R = 1500 (ref. 43). 



Figure 12. Effect of wall cooling on second mode instability in a boundary layer at 

M = 4.5 and R = 1500 (ref. 43). 
e 



Figure 13. Effect of pressure gradient on the second mode instability in a boundary 
layer at M fi = 4.5 and R = 1500 (ref. 43). 
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Figure 14. Effect of wall suction on the second mode instability in a boundary 
layer at M g = 4.5 and R = 1500 (ref. 43). 


1 / 

] /£ 


I °- 02 i v 
o.oi j y^ 


Figure 15. Variation of the maximum growth rate of secondary waves with spanwise 
wavenumber for three Mach numbers and 6 two amplitudes a of the primary 
wave at R = 1048. 8_ and an F = 83x10" for the primary wave. Dashed 
curves - a = 0.02/72 and solid curves - a = 0.01/72. Mach numbers for 
each set proceeding downward are M = 0, 0.8, and 1.6. 
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Figure 16. Variation of the growth rate of the free wave with spanwise wavenumber 
at R = 1048.8 and an F = 42.5x10“’ for three Mach numbers. 
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Figure 17. Variation of the growth rate of secondary waves with spanwise wave 
number for = 4.5; the amplitude and frequency of the primary wave are 
a = 0.01//7 and F = 120x10“ : (a) first-mode primary wave at R = 1150 

and (b) second-mode primary wave at 1950. 
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Figure 18. Variation of the growth rate of the free wave with spanwise wavenumber 
for = 4.5 and F = 60x10"°: (a) first-mode wave at R = 1150, (b) 

second-mode wave at R = 1950. 



FREQUENCY X 10«*6 


Figure 19. 


Variation of the secondary growth rate with frequency for e = 0.199, M 
= 4.5,_and R = 1950; the primary wave is a first-mode wave having an a 
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Figure 22. Pressure-coefficient distribution calculated over a backward-facing step 
whose height is -0.003 and center is at Re = 10 3 for an M = 0.5: 

IBL , — Navier-Stokes solver (grid 4). 
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Figure 23. Spatial growth rates of 2D waves having the frequnecy 80x10 

propagating over a 6 backward-facing step whose height is -0.003 


center is at Re = 10 
(grid 


for an M = 0.5: IBL, --- Navier-Stokes solver 
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Figure 24. 


Figure 25. 


Amplification factor distribution over 
height 6 is -0.003 and center is at Re = 10 
80x10" . 
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Influence of Mach number on the amplification factor of the most 
amplified 2D wave propagating over a backward-facina step whose height 
is -0.003, slope is -4.349°, and center is at Re = 10 : M ot = 0.0, 
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Figure 26. Variation of the N-factor distribution with frequency for waves 
propagating past a backward-facing step whose height is -0.003, slope is 
-4 . 349 3 . and center is at Re = 10 : M = 0.5. 




